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We report the realization of a quantum degenerate mixture of bosonic ®^Rb and fermionic ^^'^Yb atoms 
in a hybrid optical dipole trap with a tunable, species-dependent trapping potential. ®^Rh is shown to be a 
viable refrigerant for the noninteracting ^^^Yh atoms, cooling up to 2.4 x 10® Yh atoms to a temperature 
of T/Tf = 0.16(2) while simultaneously forming a ®^Rh Bose-Einstein condensate of 3.5 x 10® atoms. 
Furthermore we demonstrate our ability to independently tailor the potentials for each species, which paves the 
way for studying impurities immersed in a Bose gas. 


I. INTRODUCTION 

Mixtures of quantum degenerate gases provide a platform 
for conducting a variety of experiments, such as the creation 
of polar molecules Iffllll, the study of impurities in a Bose 
gas 0, and atom interferometry 0. In particular, there has 
been recent interest in observing the behavior of polarons on 
a lattice a and the implementation of novel lattice cooling 
schemes Em using a Bose-Einstein condensate (BEC) as a 
bath. Both of these experiments require independent tuning 
of the trapping potential for each species. While the abil¬ 
ity to create a species-dependent potential has been demon¬ 
strated in a spin-mixture of ®^Rb ll8l- [T2ll . a ^^Na-®Li mixture 
IIT 3 I and a mixture llT4ll . these methods relied on the 

light shifts from near-resonant light, which led to heating from 
large photon scattering rates. In the first case a spin-dependent 
lattice was created using the vector light shift of rubidium at 
790 nm, while in the latter two, a near resonant 670.5-nm 
(789.9-nm) beam created a large light shift for lithium (potas¬ 
sium) but a negligible one for sodium (rubidium). Mixtures 
of alkali-metal and alkaline-earth atoms GMll on the other 
hand are ideal for such experiments because their atomic spec¬ 
tra differ greatly, allowing one to use light shifts from substan¬ 
tially different wavelengths to create species-dependent traps. 

We demonstrate the realization of a degenerate mixture of 
bosonic ®^Rb and fermionic atoms. Unlike the Li-Yb 

and Rb-Sr mixtures presented in ifTSlfThll . cooling the Rb-Yb 
mixture to degeneracy has been challenging due to extremely 
unfavorable s-wave scattering lengths ifTSl between ®^Rb and 
most Yb isotopes. ^^^Yb and ^^®Yb do not form large BECs 
due to negative self-scattering lengths, and ^^^Yb phase sepa¬ 
rates from ®^Rb well before the onset of degeneracy lfT9l due 
to a large positive interspecies scattering length. While a de¬ 
generate ®^Rb-^^°Yb mixture may be possible, the s-wave in¬ 
teraction between the two species is negligible, leaving Yb 
and ^^^Yb as the only two viable candidates for a stable, in¬ 
teracting mixture with ®^Rb. 

Our mixture is trapped and evaporated in a species- 
dependent potential created by a bichromatic optical dipole 
trap (BIODT) ifTTl l20l 1^ . an ytterbium-blind magnetic 
quadrupole trap ll22ll and a rubidium-blind crossed-dipole trap 
at 423.018 nm, where the light shift for ®^Rb vanishes 1^ . 
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While this Yb isotope cannot be evaporated by itself due to 
its small s-wave scattering length of —3ao, ®^Rb serves as 
a coolant with a favorable interspecies scattering length of 
— 59ao ifTSll and we are able to sympathetically cool 2.4 x 10® 
Yb atoms down to T — O.lOTp. This represents a tenfold 
improvement in number over previous methods ll^ to cool 
^^^Yb. Eurthermore we show that this mixture is long lived 
and not significantly limited by scattering from any dipole 
beams. 

The BIODT consists of co-propagating, overlapped, 1064- 
nm (red) and 532-nm (green) dipole traps. While both beams 
provide an attractive potential for Yb, the red trap is attractive 
for Rb while the green is repulsive, allowing for independent 
control of trap depths for both species. A 423-nm crossed 
dipole beam provides longitudinal confinement for ytterbium 
atoms while the magnetic quadrupole field does the same for 
Rb atoms, giving us the ability to create species-dependent 
potentials for efficient sympathetic cooling of Yb atoms. 

Our approach, outlined in Eig. is to load laser-cooled 
nonmagnetic Yb atoms into the green dipole trap and hold 
them while a Rb magneto-optical trap (MOT) is prepared and 
transferred into the magnetic quadrupole trap for forced rf 
evaporation. The two MOTs have to be loaded sequentially 
since light-assisted collisions in a two-species MOT severely 
limit its lifetime and size. The repulsive potential created by 
the green beam prevents hot Rb atoms from heating the colder 
Yb atoms out of the dipole trap. After RE evaporation of the 
Rb to the temperature of the Yb cloud, the red beam is turned 
on, changing the repulsive green potential into an attractive 
BIODT potential and initiating thermal contact between the 
two species. Evaporation of Rb in the BIODT is then per¬ 
formed to cool the rubidium and ytterbium clouds to degener¬ 
acy. 


II. APPARATUS 

The apparatus for creating degenerate mixtures is an adap¬ 
tation of the Rb BEC system described in ll^ with a Yb 2D- 
MOT addition, similar to the setups described in EH ED. 
The two-dimensional (2D) MOT uses a field gradient of 
5.5 mT/cm (55 G/cm) provided by a pair of permanent mag¬ 
nets and is fed by an Yb oven heated to 410 °C. The 2D- 
MOT beams, detuned by —50 MHz from the ®So —tran¬ 
sition in Yb, have an intensity of 1.2/^ and are retro-reflected 
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FIG. 1. (Color online) Experimental setup for producing degener¬ 
ate mixtures of rubidium and ytterbium. The 27-^m BIODT beam 
intersects the 50-/rm wide 423-nm crossed dipole beam at an angle 
of 56°. MOT beams for both species are retro-reflected through the 
chamber windows. Gravity is along the z-axis. 


through the 2D-MOT beam viewports. A push beam detuned 
by +10 MHz from the ^Sq —>-^Pi transition accelerates atoms 
along the axis of the 2D-MOT into the experiment chamber 
where they are captured and cooled in a three-dimensional 
(3D) MOT operating on the ^Sq —line as described in 

ms). 

Using this configuration, we achieve a 3D-MOT loading 
rate of up to 6 x 10^ s“^ for The loading rate of the 

3D-MOT scales with isotope abundance for the bosonic iso¬ 
topes ^^'^Yb and ^^®Yb. However the loading 

rate for the fermions ^^^Yb and ^^^Yb is significantly worse 
than scaling with abundance due to a large differential Zeeman 
splitting between the ground and excited states Il27ll29ll . Two- 
dimensional MOTs of ^^^Yb are further complicated due to 
an unresolved hyperhne splitting in the ^Pi state. With Yb, 
we achieve a loading rate of 2.5 x 10® s-i. 

The red and green BIODT beams are sent to the experi¬ 
ment over separate optical hbers and focused at the center of 
the chamber with a waist of 27 pm, as inferred from resonant 
parametric heating measurements of Rb and ^^^Yb. Large 
mode area photonic band-gap hbers ensure a clean TEMqo 
spatial mode for each of the dipole beams and are required to 
handle the high 532-nm power typically used in the experi¬ 
ment. Conhnement along the axis of the BIODT is provided 
by the 423-nm crossed dipole trap for Yb and a quadrupole 
magnetic held for Rb. The power in each dipole beam is ac¬ 
tively stabilized. 

Atom number and temperature measurements of both 
species are performed by taking resonant absorption images 
of each atomic cloud in the y'z and xy planes shown in Fig. 

The Yb cloud is imaged on the broad ^Sq —>^Pi transition 
at 398.9 nm while the Rb cloud is imaged on the D2 line at 
780.2 nm. Dichroic mirrors placed in the imaging path sep¬ 


arate the Rb and Yb images onto different CCD cameras, al¬ 
lowing for simultaneous imaging of both clouds in both y' z 
and xy planes. 

III. EXPERIMENTAL SEQUENCE 
A. Rubidium and ytterbium MOTs 

Since ground-state Yb atoms are insensitive to magnetic 
helds, we hrst transfer Yb into the dipole trap. An 8-s MOT 
loading stage loads 2 x 10^ Yb atoms into the 3D-MOT with a 
held gradient of 2.4 G/cm. The loading rate of the Yb MOT is 
improved by spectrally broadening the 556-nm MOT beams 
from 20 kHz to 5 MHz in order to increase the MOT cap¬ 
ture velocity. A 200-ms-long cooling and compression stage 
reduces the MOT beam linewidth to 20 kHz, reduces the in¬ 
tensity from lO/g to 1.2/s, increases the held gradient 
from 2.4 to 12 G/cm resulting in a cloud with a density of 
3 X 10^° cm“^ and a temperature of 7 /iK. We then move 
the compressed MOT, using uniform magnetic helds, onto the 
focus of the BIODT before the power in the MOT beams is 
ramped down, transferring 1.5 x 10® atoms into the trap at 45 
/iK. At this stage the Yb trapping potential of the BIODT is 
provided by the 532-nm beam at 5 W and the 423-nm beam 
at 72 mW, resulting in a 325 /iK trap depth. The initial tem¬ 
perature of the Yb cloud is much larger than the 4 /tK depth 
of the crossed dipole beam, so its contribution at this stage is 
negligible. 

We load the Rb MOT from a Zeeman-slowed atomic beam 
2.5 mm above the BIODT to avoid light-assisted collisions 
heating Yb atoms out of the trap. After a 6-s MOT loading 
stage and a compression and optical pumping stage, 3.5 x 
10® atoms of Rb in the \F — l,mp = — 1) state are trans¬ 
ferred into the magnetic quadrupole trap at a field gradient 
of 192 G/cm. The held zero of the quadrupole trap is situ¬ 
ated 50 ym above the BIODT and thermalization between the 
105 /iK Rb and the 45 /iK Yb clouds is initially suppressed 
by the large repulsive potential of the green BIODT beam for 
Rb [Figjgb)]. 

B. Radiofrequency evaporation and Rb BIODT load 

The hrst stage of cooling involves forced RF evaporation 
of the Rb cloud in the magnetic trap. During this stage 
[stage (hi) in Fig. [^a)], we evaporate the Rb cloud using 
RF from 18 MHz to 5 MHz over 5 s, reducing its temperature 
to 27 /iK and increasing its density from 7 x 10^® cm“® to 
1 X 10^® cm“®. At the same time we linearly increase the 
power of the red BIODT beam to 2 W in order to gradually 
initiate thermal contact with the Yb cloud towards the end of 
this stage. Throughout this process, the Yb number and tem¬ 
perature remain constant at 1.5 x 10® and 45 /iK (Fig. [^. 

We then load the Rb cloud into the BIODT by decompress¬ 
ing the magnetic trap to 22 G/cm while further lowering the 
RF frequency from 5 MHz to 1 MHz over 2.5 s [stage (iv) in 
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FIG. 2. (Color online) (a) Experimental sequence for creating de¬ 
generate mixtures of ®^Rb and Yb from steps (i) to (v). (i) Yb 
MOT and transfer into BIODT, (ii) ®^Rb MOT and transfer into mag¬ 
netic trap (iii) RF evaporation of ®^Rb in magnetic trap, (iv) transfer 
of ®^Rb into BIODT accompanied by sympathetic cooling of Yb 
and (v) dipole evaporation of ®^Rb to degeneracy, (b)-(d) show plots 
of the full trapping potential for Rb and Yb during points and 

fs respectively. 


Fig. I^a)]. The magnetic field gradient of 22 G/cm provides 
some levitation against gravity for Rb and confinement along 
the BIODT with a trap frequency of 8 Hz. A linear ramp of 
BIODT powers, as shown in stage (iv) of Fig. increases 
the overlap between the two species and cools the Yb. Op¬ 
timal loading of Rb into the BIODT occurs when the green 


(a) 




FIG. 3. (Color online) (a) Temperature evolution of Rb and Yb 
clouds throughout the evaporation procedure. The labeled steps and 
times correspond to the ones shown in Fig|^ The evolution of the 
Yb fugacity is shown in the inset, (b) The evolution of the phase 
space density and number during the evaporation procedure. The 
steep slope of the Yb data is indicative of efficient sympathetic cool¬ 
ing by Rb. The inset shows the evolution of the Yb number over time 
during the BIODT evaporation. The total number is given by green 
circles while the number in the crossed dipole trap is represented by 
black triangles. 


power is reduced to 1.4 W and the red power to 0.8 W. We at¬ 
tribute this to a reduced capture volume at high powers, which 
stems from differing Rayleigh lengths of the 532- and 1064- 
nm BIODT beams when their waists are matched as shown in 
Fig. g The temperature of both species is reduced to 6 /iK 
at the end of decompression as shown in Fig. [^a) and the 
Yb cloud suffers minimal loss during this stage. At the end 
of stage (iv), a Rb cloud containing 7 x 10® atoms and a Yb 
cloud containing 8 x 10® atoms coexist in the BIODT at 6 ^K. 

A Ramsauer-Townsend scattering minimum at 50 fiK lIMIl 
in the interspecies thermalization cross section causes the 
^^^Yb temperature to lag behind the ®^Rb temperature during 
stage (iv). At this temperature, the interspecies cross section 
drops to 1.2 X 10® Oq, well below its low energy s-wave value 
of 4.4 X 10^ Oq. However, the thermal distribution of veloci¬ 
ties allows us to cool the Yb cloud despite the existence of a 
scattering minimum, and the two cloud temperatures converge 
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FIG. 4. (Color online) The BIODT potential for a 45 /rK deep Rb trap 
created using two different powers. The white regions are repulsive, 
(a) Rb BIODT potential for 0.8 W and 1.4 W in the 1064-nm and 
532-nm beams respectively, (b) Rb BIODT potential for 2.1 W and 
5 W in the 1064- and 532-nm beams respectively. 


to the same value towards the end of stage (iv). 


C. BIODT evaporation to degeneracy 


We perform the final stage of evaporation [stage (v) in Figs. 

and with an exponential ramp of both BIODT powers 
with a time constant r = 0.5 s over a period of 2 s. Towards 
the end of this stage the Rb cloud condenses to form a BEC 
while Yb reaches Fermi degeneracy. The Rb images show 
a strong bimodal distribution characteristic of a BEC, while 
the Yb cloud shows significant deviations from a Gaussian 
distribution (Fig. |^. After 1.7 s of evaporation in the BIODT, 
the Yb cloud reaches a temperature of T = 0.62(8 )Tf = 
220 nK with 4.2 x 10® atoms while the Rb cloud condenses 
to form a BEC of 3.5 x 10® atoms. Evaporating for 2 s cools 
the Yb cloud to T = 0.16(2 )Tf = 90 nK with 2.4 x 10® 
atoms, at the cost of reducing the Rb BEC number to 1.1 x 10®. 
This is over an order of magnitude improvement over previous 
methods to cool ^^®Yb ll24l which resulted in 8 x 10® atoms 
at r = 0.46Tp. The uncertainties in our values for T/Tp are 
dominated by noise in the number of Yb atoms and not by our 
temperature measurement. 

Near degeneracy, we extract the number and temperature of 
the ^^®Yb cloud from a two-dimensional fit of a Fermi-Dirac 
distribution. 


Li2 

n{x, y) = B — A - 


Li 2 {-z) 
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FIG. 5. (Color online) (a) A 22-ms time-of-flight (TOF) image of a 
condensed Rb cloud of 5 x 10® atoms at 112 nK, as extracted from a 
Gaussian fit to the wings of the image. The plot shows a slice of the 
image across 10 pixels, along with the fits. The bimodal fit is per¬ 
formed only on points with an OD below 3 in order to avoid artifacts 
from saturation of the probe beam, (b) A 12-ms TOF image of a de¬ 
generate Yb cloud at 90 nK with 2.4 x 10® atoms. The plot shows 
azimuthal averages of the image along with the 2D fits performed 
on the image. At this temperature the Rb cloud forms a nearly pure 
BEC, without any discernible thermal wings. 


to the absorption image, with fit parameters A, B, z, Wx, and 
Wy. Li 2 is a polylogarithm of order 2 and 2 = 
the fugacity that determines the deviation of the Fermi-Dirac 
distribution from a thermal Gaussian. For a thermal gas at 
T ^ Tp, z = 0 while in the limit T Tp the fugacity di¬ 
verges. We extract the temperature of the degenerate Fermi 
gas from Wx and Wy. The number of atoms is determined 
using two methods; through direct integration of the column 
density, and indirectly from the fugacity and temperature us¬ 


ing the relation 


Tpiio.N) ln(z)’ 

which holds in the limit T Tp. The former method 
gives a value of 2.4 x 10® atoms while the latter gives a 
value of 3.1 x 10®, using calculated Yb trap frequencies, 
(u!x,u!y,ujz) = (150,140, 75) Hz and assuming an unpolar- 
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ized Fermi gas M- 

The inset in Fig. |^a) shows the fngacity while the inset 
in Fig. I^b) shows the evolution Yb number during the final 
stage of evaporation. The fugacity rises rapidly during the last 
second of the evaporation as the cloud becomes degenerate 
but the number, which does not change significantly during 
the early stages of evaporation, drops as the trap becomes too 
shallow to support Yb against gravity. 


D. Trap lifetimes 

The lifetime of the degenerate mixture is limited by photon 
scattering of Yb atoms from the 423 nm beam with calcu¬ 
lated rate of 0.4 s“^. We measure the lifetimes of both de¬ 
generate clouds under three different conditions shown in Fig. 
1^ Rb in the absence of Yb, Yb in the absence of Rb, and a 
degenerate mixture. Scattering of 423-nm photons heats the 
Yb cloud in the absence of Rb, as shown in the inset in Fig. 
I^b). In the presence of the Rb EEC, this heating rate is re¬ 
duced and the Rb lifetime is reduced from 2.8 s to 1.5 s [Fig. 
I^a)] as the Rb cloud evaporates to keep the Yb temperature 
fixed. The Rb cloud has a peak density of 2.8(5) x 10^"^ cm“^ 
and its lifetime in the absence of Yb is limited by three- 
body recombination. Using a ®^Rb three-body rate constant 
of K 3 = 5.8 X 10“^° cm®s“^ l32l, we calculate a three body 
rate of 0.45 s“^, in approximate agreement with our measured 
lifetime. 


IV. INDEPENDENT CONTROL OE RB AND YB CLOUDS 

The flexible nature of our trap allows us to control over¬ 
lap between the two species, and independently address either 
Rb or Yb atoms. The nonmagnetic Yb atoms are unaffected 
by the magnetic trap used to confine the Rb cloud along the 
BIODT beam. Similarly the Rb cloud does not see the 423-nm 
crossed dipole beam that provides longitudinal confinement 
for the Yb gas. 

The overlap between the two clouds can be controlled by 
using a uniform magnetic field to shift the position of the mag¬ 
netic trap along the BIODT. Before BIODT evaporation, the 
magnetic Rb atoms are moved along the dipole trap while the 
position of the nonmagnetic Yb atoms is left unchanged. We 
use the Yb temperature after 0.5 s of BIODT evaporation as 
an indicator of overlap between the two species, as the non¬ 
interacting Yb atoms cannot be efficiently evaporated. This 
temperature is shown in Fig. [^as a function of the Rb cloud 
position, with a clear minimum in the Yb temperature indicat¬ 
ing optimal longitudinal overlap of the two clouds. 

Independent manipulation of the Yb cloud can be achieved 
through the 423-nm crossed dipole beam. To demonstrate this, 
we resonantly heat the Yb by modulating the intensity of the 
crossed dipole beam at the trap frequency (Fig. [^a)). This 
in turn heats up the Rb, resulting in a rise in its temperature 
as shown in Fig. [^b). When this experiment is repeated in 
the absence of Yb, we observe no heating of the Rb cloud. 



FIG. 6. (Color online) (a) Lifetime of the Rb BEC with and without 
the degenerate Yb gas present, (b) Lifetime of the Yb gas with and 
without the Rb BEC present. The inset in (b) shows the evolution of 
the Yb temperature in the presence (circles) and absence (crosses) of 
the Rb BEC. The lines in both plots are fits of an exponential decay 
to the data. 



FIG. 7. Yb temperature after BIODT evaporation, as a function 
of the position of the Rb cloud. A minimum is observed at 2.5 pK 
corresponding to the optimal longitudinal overlap between the two 
clouds. 


confirming that the two clouds are overlapped and that the 
crossed dipole beam only affects Yb. 
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FIG. 8. (Color online) (a) Yb temperature as a function of 423 nm 
modulation frequency, showing a resonance at 310 Hz. (b) Rb tem¬ 
perature as a function of modulation frequency in the presence, (red 
filled circles) and absence (black empty circles) of Yb. The Rb cloud 
is heated only in the presence of Yb confirming that both species are 
overlapped and that the 423 nm beam does not affect Rb. 


V. CONCLUSIONS AND OUTLOOK 


We have created a degenerate Bose-Fermi mixture of ®^Rb 
and Yb in a species-dependent potential. Our ability to in¬ 
dependently tune trap potentials allows us to sympathetically 
cool Yb with minimal loss, allowing us to create large degen¬ 
erate Fermi gases. We further demonstrate that the degenerate 
mixtures are stable and long lived, with the loss rate limited by 
slow photon scattering from the 423-nm crossed dipole trap. 

A species-selective optical lattice created by the 423-nm 
laser will allow us to realize several lattice cooling schemes 
Elii. Furthermore, our method of creating degenerate mix¬ 
tures is extremely flexible and can be extended to other alkali 
metals such as sodium or cesium, providing a route for a vari¬ 
ety of other degenerate alkali-metal-alkaline-earth mixtures. 
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